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Abstract 

The induced gluon radiation of a high energy quark in a finite-size QCD medium is 
studied. For a sufficiently energetic quark produced inside a medium we find the radiative 
energy loss AE q oc L 2 , where L is the distance passed by quark in the medium. It has 
a weak dependence on the initial quark energy E q . The L 2 dependence turns to L 1 as 
the quark energy decreases. Numerical calculations are performed for a cold nuclear 
matter and a hot quark-gluon plasma. For a quark incident on a nucleus we predict 
AE q « 0.1E q (L/10fm) p , with (3 close to unity. 



The radiative energy loss of a high energy parton in a QCD medium is under active 
investigation nowadays [1-5]. In classical electrodynamics the radiation of a charged 
particle in a dense medium was first considered long ago by Landau and Pomeranchuk 
|| . The quantum treatment of this phenomenon was given by Migdal J7|. In Ref. 0] 
(see also f8|) we developed a new path integral approach to the bremsstrahlung in a dense 
medium applicable both in QED and QCD. In the present paper we evaluate within the 
formalism of Ref. |§ the radiative energy loss of a fast quark, AE q , propagating through a 
finite-size uniform QCD medium. We consider both a cold nuclear matter and a hot quark- 
gluon plasma (QGP). Following [fj] we model QGP by a system of static scattering centres 
described by the Debye screened potential oc exp(— r^r,)/r, where /io is the color screening 
mass. For the screening mass we use perturbative formula /xp = (1 + n F /6y/ 2 g s T ||, 
where g s = \jAma s is the QCD coupling constant, T is the temperature of QGP. We 
assume that a fast quark produced at z = through a hard mechanism propagates in a 
medium of extent L along z axis. 

Neglecting the multigluon emission the radiative energy loss can be written as 



where E q is the initial quark energy, x is the Feynman variable for the radiated gluon, 
and dP/dx is the probability of gluon radiation as function of x. In the approach of Ref. 
0] an evaluation of dP/dx is reduced to solving a two-dimensional Schrodinger equation 
in the impact parameter space. The longitudinal coordinate z plays the role of time. This 
Schrodinger equation describes evolution of the light-cone wave function of a spurious 
three-body qqg color singlet system. The relative positions of the constituents of the 
qqg system in the impact parameter space are p q = —px, p q — 0, p — (1 — x)p. The 
corresponding Hamiltonian has the form 
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Here p(x) = E q x(l — x) is the reduced "Schrodinger mass", n(z) is the medium density, 
and o" 3 (p, x) is the cross section of interaction of the qqg system with a medium constituent 
(color centre for QGP and nucleon for nuclear matter). In the case of QGP on the rhs of 
(H) summation over triplet (quark) and octet (gluon) color states is implicit. 

In order to simplify the analysis we neglect the q — > qg spin-flip transitions which give 
a small contribution to the energy loss. Then the radiation rate is given by 0] 
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Here the generalization of the QED vertex operator of Ref. Q to QCD reads 

q g [4-4x + 2x 2 ] p(6) ■ P(6) 
9&>^x)- ^)—' (5) 

if is the Green's function for the Hamiltonian (|2]), K v is the vacuum Green's function, 
Lf = 2E q x(l — x)/[m? q x 2 + m? g (l — x)] is the so called gluon formation length (time), m q 
is the quark mass and m g is the mass of radiated gluon. The latter plays the role of an 
infrared cutoff removing contribution of the long- wave gluon excitations which cannot be 
treated perturbatively. In contrast to the expression of Ref. [Q for the bremsstrahlung 
spectrum for an electron incident on a target of Ref. in which the integration over 6 
starts from — oo, in (|J) we integrate over 6 from 6 = 0, i.e. from the point where a fast 
quark is produced by hard scattering. 

The three-body cross section entering the imaginary potential (^) can be expressed in 
terms of the dipole cross section for color singlet qq pair [|l(]], cr 2 (p), 

9 1 
<*3(p,x) = ~Mp) + a 2 ((l-x)p))} - -a 2 (xp) . (6) 

The radiation rate is dominated by the contribution from p < l/m g where 02 (p) = 
C*2(p)p 2 and C 2 (p) has a smooth (logarithmic) dependence on p JTT|, pT0f| . This allows one 
to estimate the energy loss replacing C 2 (p) by C 2 (l/m g ). Then a$(p,x) = Cs(x)p 2 , with 
C$( x ) — + (1 — x ) 2 } ~ x2 }C2({/m g )/8, and the Hamiltonian (|l|) takes the oscillator 
form with the frequency 
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Making use of the oscillator Green's function after some algebra one can represent the 
bremsstrahlung rate (II) in the form 



dP da BH 

— = Ln— — S(r],l) 

ax ax 



(7) 



where 



da 



BH 



4a s C 3 (x)(4-4x + 2x 2 



dx 9nx[m 2 x 2 + m 2 (l — x)] ' 
is the Bethe- Heftier cross section. The suppression factor S(r),l), depending on the di- 
mensionless variables 
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is given by 



S( V ,l) = SW( V J) + S^( V J), 
5« (rj, I) = A Re j dyi J dy2 exp 



sin(0y 2 ) 
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_cos(0yi) (taa(^yi) + 4>y 2 ) 

with = = exp(— in/ 4). The two terms on the rhs of (pTT| ) correspond in (fj) to the 

contributions from the integration regions £i < £2 < L and £1 < L < £ 2) respectively. The 
variables in (|T2"D, ( |T3"D in terms of those in (|J) are yi = (L — £i)|fi|, 2/2 — (£2 — £i)|^| ( m 
(0)) and 1/2 = (£2 — L)\Q\ (in (pT3|)). In arriving at ( [[3D we have used representation of the 
first Green's function in the square brackets in @ through convolution of the oscillator 
Green's function (for the interval (£i,L)) and the vacuum one (for the interval (L, £2))- 
Notice that the functional form of our results at x <C 1 differs from the one obtained in 
H within the soft gluon approximation. 



In a medium it is either Lf or l/|f2| which sets the effective medium- modified formation 
length L'j = min(Lj, 1/|0|), which is the typical value of £2 — £1 in (|j) for L ^> L'j. The 
finite-size effects come into play only at L < Ly, i.e. I <h = min(l, 1/rj). ^From ([11 



(|13|) we find S(r],l) ~ — Z 2 logZ as I — > 0. The source of this suppression of radiation at 
small L is obvious: the energetic quark produced through a hard mechanism loses soft 
component of its gluon cloud and radiation at distances shorter than the time required 
for regeneration of the quark gluon field turns out to be suppressed. For / ^> Z S(r], I) 
reduces to that for the infinite medium, for which S(r/,l = 00) w 3/r]^/2 ( 77 ^> 1 ) and 
S(r), I = 00) « 1 — 16?y 4 /21 (77 <C 1) were derived in [Q. Notice, that according to (Q), 
(|10|) 7] —>■ and / — >• 00 as a; — > 0, 1 and the Bethe-Heitler regime takes place in these 
limits. 

Before presenting the numerical result, let us consider the energy loss at a qualitative 
level. We begin with the case of a sufficiently large E q such that the maximum value of 
Lj, L'j(max), is much bigger than L. Taking into account the finite-size suppression of 
radiation at L'j > L, we find that AE q is dominated by the contribution from two narrow 
regions of x: x < 5 g ~ Lm^/2l E q and 1 — x < 5 q ~ Lm^/2l E q . In both the regions 
the finite-size effects are marginal and the energy loss can be estimated using the infinite 
medium suppression factor. For instance, 

16a s C 3 (0)E q Ln 
2 
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Using (H) one can show that r](x < 5 g ) < 1 at L < m 2 g /2nCz{Q). In this region of L in 
(|14]) we can put S(r](x),l = 00) ps 1 and find AE q ~ 0.25a s C 3 (0)nL 2 , which does not 
depend on the quark energy. At L ^> m 2 g /2nC^{Q) the typical values of rj in (0) are 
much bigger than unity, and using the asymptotic formula for the suppression factor we 
obtain AE q ~ a^Cs^nL 2 . Similar analysis for x close to unity gives the contribution 
to AE q suppressed by the factor ~ 1/4 as compared to that for small x. Notice that in 
this L 2 regime, despite the l/m 2 q infrared divergence of the Bethe-Heitler cross section, 
AE q has only a smooth m 9 -dependence originating from the factor C3. We emphasize 



AE q (x <S g ) ~ l&asC ^ qLn JdxS( V (x),l = 00) . (14) 



that the above analysis of the origin of the leading contributions makes it evident that L 2 
dependence of AE q cannot be regarded as a consequence of the Landau-Pomeranchuk- 
Migdal suppression of the radiation rate due to small angle multiple scatterings. 

The finite-size effects can be neglected and AE q becomes proportional to L if 
LJmax) <C L. If in addition the typical values of 77 are much bigger than unity, from 
(E)> (0)> (H) along with the asymptotic form of S(r], I = 00) at 77 ^> 1 one can obtain the 



following infrared stable result AE q rs 1.1ol s LJ nCz(<S)E q . 



In numerical calculations we take m g = 0.75 GeV. This value of m g was obtained in 



1 12] from the analysis of HERA data on structure function F2 within the dipole approach 
I3| to the BFKL equation. It is also consistent with the nonperturbative estimate [|14[ 



of the gluon correlation radius in QCD vacuum. For scattering of the qqg system on a 



nucleon, we find from the double gluon model [|TTJ C2(l/m g ) ~ 1.3 — 4 where the lower 
and upper bounds correspond to the t-channel gluon propagators with mass 0.75 and 
0.2 GeV, respectively. The latter choice allows one to reproduce the dipole cross section 
extracted from the data on vector meson electroproduction [[0J. However, there is every 
indication T3| that a considerable part of the dipole cross section obtained in [|TjJ 
comes from the nonperturbative effects for which our approach is not justified. For this 
reason we take C2(l/m 9 ) = 2 which seems to be plausible estimate for the perturbative 



component of the dipole cross section |T2|]. For scattering of the qqg system on color 
centre we estimated C*2(l/m 9 ) using the double gluon formula with the Debye screened 
gluon exchanges. For T = 250 MeV we obtained C2(l/m 3 ) ~ 0.5 for triplet centre. For 
octet centre the result is Ca/Cf = 9/4 times larger, here Ca{Cf) is the octet (triplet) 
second-order Casimir invariant. For quark mass, controlling the transverse size of the qqg 
system at x ~ 1, we take m q = 0.2 GeV. Notice that our prediction for AE q is insensitive 
to the value of m q . 

For nuclear matter taking n = 0.15 fm~ 3 and a s = 1/2 for L < 5 fm we obtained 
AE q w a(L/5fm) /3 , with a w 0.55, 1, 1.23 GeV and (3 w 1.5, 1.85, 1.95 for E q =10, 50, 
250 GeV. Calculations with a s = 1/3 for QGP at T = 250 MeV yield for the same energies: 



5 



a ~ 4.2, 10.2, 14.8 GeV and (3 ~ 1.2, 1.65, 1.9. The above values of (3 were determined for 
L < 5 fm. In the region 5 < L < 10 fm (3 is by 10-20 % smaller. At E q > 250 GeV a and (3 
flatten. Notice that L'^rnax) ~ 5 — 10 fm for E q ~ 10 — 40 GeV in nuclear matter, and for 
E q ~ 150 — 600 GeV in QGP. Thus our numerical results say that the onset of L 2 regime 
takes place when L'j(max) / L > 2. The closeness of (3 to unity at E q = 10 GeV for QGP 
agrees with a small value of L'f(max) (~ 1 fm). We checked that variation of m q gives a 
small effect. The m 9 -dependence of AE q becomes weak at E q > 50 GeV. However, it is 
sizeable for E q ~ 10 - 20 GeV. For instance, AE q {m g = 0.375)/ AE q (m g = 0.75) ~ 1.5 at 
E q = 10 GeV, L ~ 5 fm. Our predictions for AE q must be regarded as rough estimates 
with uncertainties of at least a factor of 2 in either direction. Nonetheless rather large 
values of AE q obtained for QGP indicate that the jet quenching may be an important 
potential probe for formation of the deconfmement phase in AA collisions. 

We also studied the energy loss of a fast quark incident on a target. In this case the 
radiation by initial quark is allowed and the lower limit of integration over £i in (|4]) must 
be replaced by — oo. For the bremsstrahlung in QED this situation was discussed in ||. It 
was shown that after expanding the medium Green's function in a series in the potential 
the spectrum can be represented as a sum of the Bethe-Heitler term and an absorptive 
correction. For our choice of the gluon mass the absorptive correction is relatively small. 
This means that AE q oc E q Lna s Cs(0) j 'm 2 ,. For nuclear matter in the region L < 10 fm 
the numerical calculations give AE q m 0.1E q (L/10 fm)^ with [3 m 0.9 — I for E q < 50 
GeV and (3 ~ 0.85 — 0.9 for E q > 200 GeV. This result differs drastically from prediction 
of Ref. m AE q pa 0.25(L/lfm) GeV. Our estimate is in a qualitative agreement with the 
longitudinal energy flow measured in hard pA collisions with dijet final state |16| and the 
energy loss obtained from the analysis of the inclusive hadron spectra in hA interactions 

I would like to thank D. Schiff for discussions and hospitality at LPTHE, Orsay, where 
this work was completed. 
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